Harmonically trapped dipolar fermions in a two-dimensional square lattice 
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We consider dipolar fermions in a two-dimensional square lattice and a harmonic trapping po- 
tential. The anisotropy of the dipolar interaction combined with the lattice leads to transitions 
between phases with density order of different symmetries. We show that the attractive part of the 
dipolar interaction results in a superfluid phase which is suppressed by density order. The trapping 
potential is demonstrated to make the different phases co-exist, forming ring and island structures. 
The phases with density and superfluid order can overlap forming regions with supersolid order. 



The trapping of dipolar atoms and molecules is a 
promising new research field. The anisotropy of the 
dipole interaction offers unique opportunities for explor- 
ing novel few-body [TJ [2] and many-body quantum sys- 
tems [21 H]. Experimentally, one has realized Bosc- 
Einstein condensates of 52 Cr atoms [SJ E] and 164 Dy 
atoms |7] with large magnetic dipole moments, as well as 
gases close to quantum degeneracy of 40 K 87 Rb molecules 
with an electric dipole moment [5]. Furthermore, the first 
experimental steps toward realizing dipolar molecules in 
an optical lattice have recently been reported ITU] . 
The lattice makes the physics very rich: Density ordered 
phases with a complicated unit cell pj], liquid crystal 
phases [12], and a supersolid phase pR|] have been pre- 
dicted to exist for fermionic dipoles in a 2D lattice with 
dipole moments perpendicular to the lattice plane, and 
tilting the dipoles towards the lattice plane leads to bond- 
solid order and p-wave superfluidity at half-filling |14) . 

We consider fermionic dipoles in a 2D square lattice at 
zero temperature. A harmonic potential, which is always 
present in trapped atomic/molecular systems, is included 
exactly, since the characteristic lengths of the ordered 
phases can be comparable to the system size for experi- 
mentally realistic systems. This means that one cannot 
simply resort to the local density approximation. A main 
purpose of the present paper is to study the rich physics 
coming from the interplay between the anisotropic dipole 
interaction, the optical lattice, and the inhomogeneity in- 
duced by the trapping potential. For experimentally re- 
alistic systems, we demonstrate the existence of compet- 
ing phases with density order of checker-board or stripe 
symmetry, and superfluid order. Due to the trapping 
potential, these phases can co-exist and sometimes even 
spatially overlap leading to regions with supersolid order. 

Basic formalism The fermions have a mass m and a 
dipole moment d which is aligned by an external field 
forming the angle Op with the z-axis perpendicular to 
the lattice (xy) plane and the azimuthal angle <pp with 
respect to one of the lattice vectors parallel to the x- 
axis, see Fig. [l] In the lowest band approximation, this 
system is described by the extended Hubbard model with 



the Hamiltonian H — Hq + V where 

(ij) i V J 

and 

V =\Y. V ^)^ ( 2 ) 

with Tij = Vi —Yj. Here, removes a dipole at site i with 
position Yi, fii — c\ci, fi is the chemical potential, t is the 
hopping matrix element between nearest neighbor sites 
(ij), and i ^ j since we are considering identical fermions. 
The trapping frequency is u> and the interaction between 
two dipoles separated by r is 




FIG. 1: (color on-line) We consider dipoles in a square 2D 
lattice in the a:t/-plane. The dipoles are aligned forming an 
angle 8p with the z-axis and the azimuthal angle <f>p with the 
z-axis which is parallel to one of the lattice vectors. 
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Since the dipolar interaction ([3]) has both attractive 
and repulsive regions, the system exhibits both pair- 
ing and density instabilities depending on (9p,<f)p). To 
model this complex behavior, we decouple the interaction 
V using the mean-field approximation which yields 




where (cyCj) is the pairing order parameter. Even though 
fluctuations are important in 2D, we expect mean-field 
theory to capture the existence and competition between 
different ordered phases at T = 0. Indeed, mean-field 
theory is widely used in the high-T c community to de- 
scribe the competition between e.g. anti-ferromagnetic 
and superffuid ordering 15 . We have not included the 
Fock term in ([4]), since one dipole interacts with many 
others making the problem similar to a high dimensional 
one, for which the Hartrec term dominates [16 . The 
Fock term has recently been shown to lead to bond-solid 
order phases at half-filling [TJ] . 

The mean-field Hamiltonian H + Vmf is diagonalized 
by the Bogoliubov transformation Cj = > o( u r)% + 
v r]*lri) where 7^ are fermionic operators annihilating a 
quasi-particle with energy E n . The wave functions ut 
and satisfy the Bogoliubov-de Gennes equations 

-MS)- (5 > 

with Aij — VD{^ij){cjCi) and 
Lij = -tSfa) + V D (r lk )(n k ) + ^-^ 2 r? - (6) 

k 

Here Sij and 6^ are the Kronecker delta functions con- 
necting on-site and nearest neighbor sites, respectively. 
Self-consistency is obtained iteratively through the usual 

relations: (hi) = J2e >q Kl 2 and i^c,) = J2e v >o u I, v v* ■ 
No trapping potential Consider first a system with no 
trapping potential. We plot in Fig. [2] the ground state 
as a function of the dipolar angles (9p,(j>p) in the limit 
of strong interactions g/t — 16 with g = D 2 /a? where 
a is the lattice constant. The filling fraction is / — 
iV _1 ^,-(ni) = 1/2 with N the number of lattice sites. 
We consider two possible ground states: One with a 
checker-board density order and one with a striped den- 
sity order. The + symbols and the x symbols indicate 
when mean-field theory on a 26 x 26 lattice predicts the 
checker-board phase and the striped phase to have the 
lowest energy, respectively. In order to take into account 
the long-range nature of the dipolar interaction we have 
calculated the Hartree term by duplicating the 26 x 26 
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FIG. 2: (color on-line) Top: The phase-diagram for an un- 
trapped system for strong coupling g/t — 16 and half-filling. 
The dashed line shows the boundary between the checker- 
board and the striped phase. Bottom: The nearest neighbor 
pairing \(c i+£ Ci)\ (x's) and stripe order \(hi — n i+y )\ (+'s) 
as a function of g/t with (6p,(f>p) — (n/2,0). The o's show 
\{rii — ni+y)\ when we take the gap to be zero. The V's show 
\{&i+xCi)\ when the density is taken to be homogeneous. The 
inset shows a cut of (c i+ xCi) along the x-axis when g/t — 0.9. 



lattice so that it constitutes a lattice of 9 • 26 x 9 • 26 
sites. For small 9p, the checker-board phase is favored 
whereas for larger 9p & phase with stripes along the IE- 
direction has the lowest energy. This is easily understood 
for <fip = 0, where the dipoles are aligned head-to-tail in 
this phase which clearly minimizes the interaction energy. 
When (j>p > 0, the alignment is not perfect by stripe for- 
mation along the x-axis, but the interaction energy is 
still minimized although it requires larger tilting angles 
as expected. This is illustrated further by the dashed 
line separating the checker-board phase from the striped 
phase, which is a result of comparing the classical inter- 
action energy in the two phases. The good agreement 
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between the numerics and this calculation demonstrates 
that for strong coupling g/t ^> 1 the kinetic energy can 
be neglected, and the problem becomes classical. 

When 8p > arcsin(l/v / 3) ~ 0.27T, the interaction has 
attractive regions, and we now examine the competi- 
tion between the resulting pairing instability and the 
instability towards density order. In Fig. [2j we plot 
as x's the largest value of the nearest neighbor pairing 
|(ci+£Ci}| as a function of the coupling strength g/t with 
x a unit vector along the x-direction. We have chosen 
(9p,(f>p) = (ir/2,0) and / = 1/3. The calculations are 
performed on a 27 x 27 lattice. For weak coupling, the 
ground state is a superfluid. The cross section along the 
x-axis of the pair wave function (cjc/) with = (0,0) 
plotted in the inset, illustrates that it is odd under inver- 
sion as for the homogeneous case [171 HH] . For simplicity, 
we refer to this as p-wave symmetry in the following, even 
though the pair wave function in general contains higher 
odd components of angular momenta. For stronger cou- 
pling, the pairing vanishes as the striped phase emerges. 
The stripe order defined as (fii — n i+ y), with y a unit vec- 
tor along the y-direction so that {h/) and (fii+y) give the 
densities for a site in the stripe and next to the stripe re- 
spectively, is plotted as +'s in Fig. [2} We also plot as o's 
the stripe order parameter when pairing is not included 
in the calculation, and as V 's the pairing order parameter 
when the stripe formation is not included. This demon- 
strate that the stripe order is insensitive to pairing since 
the critical coupling strength for the formation of stripes 
essentially does not change when pairing is included. The 
pairing on the other hand does not vanish with increas- 
ing coupling if the stripes are suppressed. We therefore 
conclude that it is the pairing which is suppressed by the 
stripe formation and not the other way around. 

More complicated density order with larger unit cells 
can appear for untrapped systems, as has recently been 
demonstrated for Op = [IT]. Here, we focus on the ex- 
perimentally most relevant orders with the smallest unit 
cell, since the trapping potential will complicate the ob- 
servation of orders with a large unit cells. 
Trapped case We now examine the interplay between 
the harmonic trapping potential and the competition be- 
tween density and pairing instabilities. 

Angle 9p = Consider first the case when the dipo- 
lar orientation is perpendicular to the lattice so that the 
interaction is purely repulsive. Figure [3] shows the den- 
sity profile for trapped dipoles with Cj = ujay^ m/t = 0.24 
which is a relatively weak trapping potential so that there 
are regions with phases resembling those for the case with 
no trapping potential. For g/t = 0.5 and 144 particles, 
there is no density order whereas for stronger coupling 
g/t = 1 and 207 particles, the system exhibits a checker- 
board density profile in the center of the trap. In this 
region, / ~ 1/2 which is optimal for the checker-board 
phase |llj . The checker-board phase is surrounded by a 
normal phase with a lower density. Experimentally, sim- 
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FIG. 3: (color on-line) Left: The density profile through the 
center of the trap. The blue solid line is for g/t — 1 and 
207 particles trapped and the green crosses are for g/t — 0.5 
and 144 particles trapped. The red dotted line indicates the 
trapping potential in units of 2.5t with Q — 0.24. Right: 2D 
plot of the density in the lattice plane for the g/t = 1 case. 
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FIG. 4: (color on-line) Left: The particle density for 205 
dipoles with 6p = tt/2, g/t = 0.85, and u = 0.11. Right: The 
nearest neighbor pairing (\{c i+ii Ci)\ + \(c i - i c i )\)/2. 

ilar shell structures have been observed for atoms with a 
short range interaction [T5] . 

Angles (9p,(f)p) — (tt/2, 0) Consider next the case 
when the dipoles are aligned in the plane along the x- 
axis (<fip = 0). There is then a competition between 
density and pairing order. In Fig. [4j we plot the density 
and the nearest neighbor pairing ( | (ci+ x Ci) | + | (ci-xc/) |)/2 
(symmetrized to reduce trap effects) for g/t = 0.85, w — 
0.11 and 205 dipoles trapped in a 39 x 39 lattice. For 
this set of parameters, there is no density order since 
the density is low and the interaction weak. This results 
in p-wave pairing throughout most of the cloud. The 
cloud profile is slightly elongated in the x-direction due 
to the anisotropy of the dipolar interaction in analogy 
with what has been observed for dipolar condensates [1] . 

Figure [5] depicts the density and the nearest neighbor 
pairing for a stronger coupling strength g/t = 1 with 
w = 0.11 and 180 dipoles trapped on a 39 x 39 lat- 
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FIG. 5: (color on-line) Left: The density profile for 180 dipoles 
with 6p — 7r/2, g/t = I, and Co — 0.11. Right: The nearest 
neighbor pairing (\{c l+i Ci) \ + |(c 1 _ 4 c i )|)/2. 



tice. There is a pronounced stripe order in the center 
of the trap with an average filling fraction / ~ 1/2 which 
squeezes the pairing away from the center into two is- 
lands centered at x = 0. This intriguing island struc- 
ture is a consequence of the anisotropy of the interac- 
tion. Since the interaction is attractive when the dipoles 
are aligned head-to-tail (x-direction) and repulsive when 
they are side- by-side (y-direction) , the average interac- 
tion for a given radius is attractive in the regions close to 
the j/-axis whereas it is repulsive in the regions close to 
the x-axis. Thus, pairing can exist in the islands around 
x = away from the center where stripe order dominates, 
whereas it is suppressed in the regions around y = 0. 
These islands of pairing should be compared with the 
ring structures predicted for a two-component trapped 
fermi gas with a short range isotropic repulsive inter- 
action where anti-ferromagnetic order competes with d- 
wave pairing [20] . Remarkably, the stripe order is not 
completely suppressed in the two islands of pairing. The 
pairing in fact oscillates in-phase with the stripe order. 
This indicates that the trapping potential induces regions 
of pairing co-existing with density order as in a supersolid 
- a phase which has been subject to intense investigations 
since its theoretical prediction long ago [2lT - f25] . 

The lattices considered here have experimentally real- 
istic sizes, and with electric dipole moments up to several 
Debye for the experimentally relevant KRb, RbCs, and 
LiCs molecules |5I [2l)U2"5] . one can easily reach the strong- 
coupling regime with g/t ^ 1 using typical values for an 
optical lattice. The density ordering can be observed di- 
rectly by in-situ imaging |29] or by time-of-flight experi- 
ments which also can detect pairing correlations EU] • 

In conclusion, we have shown that dipolar fermions at 
T = in a 2D square lattice exhibit phases of density 
order with different symmetries as well as a p-wave su- 
perfluid phase. The system is unstable towards pairing 
when the interaction has attractive channels. However, 
any superfluidity is suppressed if there is an instability 



towards stripe formation such as in regions with a den- 
sity close to half filling. The trapping potential leads 
to an inhomogeneous filling fraction which results in the 
co-existence of several phases. These phases can overlap 
resulting in the presence of regions with supersolid order. 
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